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We demonstrate storage and manipulation of one qubit encoded into a decoherence-free subspace (DFS) 
of two nuclear spins using liquid state nuclear magnetic resonance (NMR) techniques. The DFS is spanned 
by states that are unaffected by arbitrary collective phase noise. Encoding and decoding procedures reversibly 
map an arbitrary qubit state from a single data spin to the DFS and back. The implementation demonstrates 
the robustness of the DFS memory against engineered dephasing with arbitrary strength as well as a substantial 
increase in the amount of quantum information retained, relative to an un-encoded qubit, under both engineered 
and natural noise processes. In addition, a universal set of logical manipulations over the encoded qubit is also 
realized. Although intrinsic limitations prevent maintaining full noise tolerance during quantum gates, we show 
how the use of dynamical control methods at the encoded level can ensure that computation is protected with 
finite distance. We demonstrate noise-tolerant control over a DFS qubit in the presence of engineered phase 
noise significantly stronger than observed from natural noise sources. 



I. INTRODUCTION 



The ability to effectively protect the coherence properties of a quantum information processing (QIP) device against the 
detrimental effects of environmental interactions is a prerequisite for realizing any potential gain of quantum computation and 
quantum information theory [|l]]. Approaches based on noiseless (or "decoherence-free" |^]) coding offer a promising venue for 
meeting the challenge of noise-tolerant QIP. The theory of decoherence-free subspaces (DFSs) has been the focus of intensive 
development particularly by Zanardi, Lidar, and coworkers [ 3|; ?|| . Recently, the DFS idea has been incorporated within the more 
general approach based on noiseless subsystems (NSs) 13 1, which recover DFSs and their benefits as special instances. 

The primary motivation behind "passive" noise control strategies relying on either DFSs or NSs is to take advantage of specific 
symmetries occuiTing in the noise process to single out subspaces or subsystems of the physical information processor that are 
inaccessible to noise. Once information is appropriately encoded into such noiseless structures, robust storage is ensured without 
requiring further active coiTection - as long as the underlying symmetry dominates. These features, together with their stability 
against symmetry-perturbing errors [^,||j8| and the consequent potential for concatenation with quantum error-coiTecting codes 
10], make noiseless codes natural candidates as robust quantum memories. To date, experimental implementations include 
studies of DF states in quantum optical systems [|T^] , and one-bit quantum memories based on both a DFS of two trapped ions 
Jl^ ] and a NS of three nuclear spins JT^ ] . 

Achieving robust quantum information storage represents only a first, though indispensable, step toward the goal of reliable 
QIP. An important advance in this direction came from the identification of universality schemes, which in principle enable 
DFSs (or NSs) to support universal encoded quantum computation in a way that remains fully protected against noise. Both 
existential [[Tl 17 1 and constructive results Jig ] have been established. While the latter are especially appealing for a class 
of proposed quantum computing architectures governed by Heisenberg exchange interactions [|l9[|, implementations of these 
schemes remain difficult due to the stringent symmetry and tunability requirements on the control Hamiltonians. 

Here, we take a first experimental step towards encoded quantum computation by demonstrating universal control over a 
one-bit DF quantum register of two nuclear spins. A novel key ingredient we use to irnplement encoded quantum gates is the 
combination of robust control design with the use of dynamical decoupling methods po|-p3|] directly on encoded degrees of 
freedom. Our results suggest that this may serve as a useful strategy for practically coping with the constraints required for DF 
computation. 

The paper is organized as follows. In Sect. II we review the collective decoherence model that is relevant to the work, along 
with the prescriptions from the DFS theory for both protected storage and manipulation of quantum information in a two-qubit 
system. In Sect. Ill, we outline our proposed approach to noise-tolerant control of DFS encoded qubits based on concatenating 
encoded decoupling methods with robust control design. The general principles are developed starting from the physical NMR 
setting relevant to the experiment. Sect. IV contains an account of the control techniques used in the experiment and the 
reliability measures adopted to quantify the accuracy of the implementation. In particular, a notion of gate entanglement fidelity, 
generalizing Schumacher's definition to allow a desired unitary evolution on the quantum data, is proposed and related to other 
fidelity metrics relevant to QIP. The experimental results demonstrating protected storage and universal protected quantum logic 
are presented and discussed in Sect. V and VI, respectively. 
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II. PROTECTING QUANTUM INFORMATION AGAINST COLLECTIVE DECOHERENCE 



A. Collective decoherence 



For a system S composed of n qubits, a purely decohering, collective interaction arises when the qubits couple symmetrically 
to a single environment E and no exchange of energy takes place between S and E. Physically, this model accounts for relaxation 
due to fully correlated fluctuations of the energy levels of each qubit - a situation that is approached if the qubits are close enough 
relative to the correlation length of the environmental coupling and the latter commutes with the natural Hamiltonian. Although 
not always applicable, this decoherence model has a practical significance for QIP. In particular, collective dephasing was shown 



to play a major role in quantum devices based on trapped ions [15|. In NMR systems, dephasing caused by fully correlated 
fluctuations of the local magnetic field provides the dominant relaxation mechanism of quantum coherences between identical 
species in sufficiently small, rigid molecules [^4[|. 

If _ff = Hs O 1 _E + Is® He + Hs e represents the Hamiltonian for the joint system plus environment, collective phase 
damping corresponds to an interaction Hamiltonian of the form 

HsE = Jz®B,, (1) 

where the operator Bz acts only on E, and measures (in units h/2) the projection of the total spin angular momentum along 
the quantizing axis z. Thus, Jz — '^z terms of the Pauli operator acting on the jth qubit. Starting from an initial state 

Qin = Y^in) (^iril, the systcm alone evolves after a time t into Qin ^ Qout ~ S{gin), where S is the super-operator associated 
with H nin. Under the assumption that S and E are initially uncorrelated, Qout can be expressed as 



Qout 



for a set of Kraus operators {Fa} [|25|]. Because [Hs, J^] = for pure decoherence, the unitary contribution due to Hs can be 
separated out as Fa = Ea exp{—iHst), leaving a set of error operators {Ea}, that still satisfy J^a FlEa = is, and describe 
the non-unitary effects of the environment. If t is sufficiently short (or, equivalently, the coupling is weak enough), the operators 
Ea can be expressed as linear combinations of the basic error generator Jz and the identity. For arbitrary collective decoherence, 
the possible errors that the coupling ([l|) can induce belong to the interaction algebra [[l^] Az generated by Jz i.e., the algebra 
containing all the linear combinations of arbitrary powers of Jz and I5. By construction, Az is Abelian, expressing the fact that 
pure decoherence is energy-conserving - thus diagonal in the computational basis. 

We focus on n = 2 qubits, in which case every element in Az can be written as a linear combination of three operators, 
t,Jz =cr]+cr2^ J2- equivalently, we replace the latter with 0-1 cr2 = J^/2-lil = I5 henceforth). If {|00), |01), |10), |11)} 
denotes the computational basis, an equivalent choice as a basis for Az are the three orthogonal projectors on subspaces with 
definite z-angulai- momentum, n+2 = |00)(00| = (1 + + cr]cr2)/4, Hq = |01)(01| + |10)(10| = (1 - (tI<tI)/2, n_2 = 
|11) (11| ~ (1 — + (tJ(t^)/4. In addition to labelling the basis states, the Jz q uantum number also allows classification of the 
possible transitions between these levels, via the so-called coherence order [ |24( | . For a quantum coherence between two states 
with i-angular momentum k, £, corresponding to the off-diagonal density matrix element |fc) {£\ (k, £ being now measured in 
units h), the coherence order mki = |A: — ^|. The behavior of the zero-, one-, and two-quantum coherence orders present in a 
two-spin system undergoing collective decoherence can be described by obtaining an explicit set of error operators {Ea}. This 
is done starting from a unitary representation Use for the joint SE evolution [p5[], 

|00)|e)^ |00)|eo), 
|01)|e)^|01)|ei), 
|10)|e)^|10)|e;), 
|ll)|e)->|ll)|e2), 

where the |e), \ek) are generally non-orthogonal environment states, and the collective nature of the interaction sets {ei\e'i) = 1. 
The overlaps (eo|ei) = (ei |e2) — e~'' , (eole^) = e~'' parametrizes the decay rates for single and double quantum coherences, 
respectively. By letting 7' = 47 (7 > 0) and by using an orthonormal basis \^a) obtained from the jca), a = 0, 1, 2, a 
choice of Kraus operators is given byi?a = {^ia\UsE\e) i.e.. 

En = n+2 + e-'Ho + e-47n_2 , 



El = VI - e-^T Up + e-T( l + e-^^)Vl -~ e'^T n_2 , (2) 

E2 = (1 - e^27)VTT^=2Tn_2 . 

An equivalent representation was derived for dephasing quantum operations engineered via gradient-diffusion techniques in 
NMR [M. In the limit of arbitrarily strong (or "crusher") interaction, 7 ^ cx), Eqs. (0) simplify to Eq = n+2, Ei = XIq, i?2 = 
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n_2, leading to the full suppression of single and double quantum coherences. However, regardless how strong, collective 
decoherence produces no decay of the zero-quantum subspace spanned by { |01) , 1 10) }. 

Zero-quantum coherences and their properties have been long appreciated in NMR, with important applications in both high- 
resolution spectroscopy in inhomogeneous magnetic fields and contrast enhancement in magnetic imaging [ ^^ . Within NMR 
QIP, zero-quantum coherences are revisited in view of their natural potential to encode protected quantum information. 



B. Decoherence-free encodings 

In the DFS approach, the first step is ensuring that the quantum data to be protected is encoded into a DFS. Mathematically, a 
DFS is a subspace of the system's state space spanned by a set of degenerate eigenvectors of all the error generators appearing 
in HsE- For global dephasing on n qubits as in dll), let H(J-) be the eigenspace corresponding to the eigenvalue jz of J^, 
jz = 71, 71 — 2, ... , — ri + 2, — n. Then states in Ti^^" remain invariant under the environmental coupling, 

J.I^l) VlV^i) e , (3) 

which also implies a degenerate action of each error operator on the subspace: 

Eali^L) = fa\^^L) , y\l^L) & U'^'^K ^Ea^Az, (4) 

for some coefficients fa — fa{jz) fulfilling |/ap — 1. If the system is initialized in a state gin = G Ti'^^'K then 

a a 

i.e., the evolution remains unitary within each TiS^'^ Thus, each TiS^'^ is a DFS under collective decoherence. The amount of 
quantum information that a given DFS is able to protect is determined by its dimension rij^ - which is simply the degeneracy 
of the corresponding -eigenvalue. In particular, for n even, the largest DFS is supported by the zero-quantum subspace 7i'^°), 
with no = n\/{n/2\Y. 

For n — 2 spins, 1-6^^ is doubly degenerate, hence it provides the smallest DFS capable to protect one qubit against collective 
decoherence. The robustness of this two-spin zero-quantum subspace under Az has been explicitly derived above. In terms of 
basis states, our DFS qubit is defined by the encoding 

co|0l)+ci|1l)=co|01) + ci|10), T^l = Ti^") = span{|OL), |1l)} , (5) 

for arbitrary complex coefficients cq , ci . It is worth stressing that the existence of a DFS is tied, at the physical level, to the 
occurrence of symmetries in the noise process. The way the latter reflect into the state space of the system both determines 
the possibility of invariant states as in (^-(Q) and the associated degeneracies. For any interaction which is diagonal in the 
computational basis, the underlying "axial" symmetry ensures that the individual a'l are conserved quantum numbers. However, 
it is only under the additional permutation symmetry characterizing collective interactions that the degenerate conserved quantum 
number Jz arises - signaling the presence of a protected structure. The emergence of degenerate degrees of freedom preserved 
under the noise remains the key ingredient for more general DFSs [p|,p|p^ and, in still more elaborated forms, NSs as well 
[ jl0| , [l2| , p^ . Taking advantage of the existing symmetries translates into major gains toward achieving noise-protected QIR For 
instance, the simple 2-bit encoding (Q) preserves a qubit against collective dephasing of arbitrary strength, to be contrasted with 
independent phase errors - where protection can be achieved only with finite distance using a quantum error-correcting code. 



C. Decoherence-free manipulations 

Once protected storage of quantum information is obtained, the next step is to ensure that universal quantum gates are im- 
plemented without ever leaving the DFS. Because the states spanning a DFS are characterized by precise symmetry properties, 
symmetries are likewise crucial in determining the control operations to be applied for effecting DF quantum logic. Clearly, the 
allowed gates must map DFS states to DFS states. However, as any physical gate takes a finite time to execute, invoking unitary 
manipulations that preserve the DFS at the conclusion of the gate is not sufficient. To guarantee that the system remains within 
the DFS during the entire gating time requires the stronger condition that gates are generated by Hamiltonians that themselves 
respect the symmetry. Thus, the general problem requires identifying a universal set of Hamiltonians which satisfy the correct 
symmetry constraints and involve at most two-body interactions [|l8|]. 

In our case, because a single DFS qubit is involved, this universal set of control Hamiltonians is composed of two observables 
generating an encoded u(2) Lie algebra, exponentiation then giving the whole group U(2) of encoded one-qubit transformations. 
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A necessary and sufficient condition for an Hamiltonian A to preserve the zero-quantum DFS can be obtained by demanding 
that (00|yl(co|01) +ci|10)) = 0, (ll|^(co|01) +ci|10)) = for arbitrary cq, ci. This leads to the following matrix form for A 
with respect to the computational basis: 



/ ai c 

02 6 

5* 03 

\c* 04 



(6) 



for real coefficients , j = 1 , . . . , 4, and possiblycomplex b,c. In particular, this includes all the hermitian operators belonging 
to the so-called commutant of the error algebra [ [lO|jl^ ], A'^ — {X : [X, Jz] ~ 0}, which collects all operators commuting with 
the noise. Because every operator in A'^ can be represented as a linear combination of the identity, the one-bit operators a^, and 
the two-bit couplings ct^ct^, ■ — ct^ct^ + cr^cr^ + o'lo'l (Heisenberg coupling), Hamiltonians in A'^ have c = - implying 
that all the DFSs are in fact preserved. An additional constraint (so-called independence [p^) can be imposed on Hamiltonians 
in A'z by also requiring oi = 04 = 0, in which case A has zero entries outside the selected zero-quantum DFS. 

With respect to the DF encoding (^, a choice of operators that act as independent, encoded tr^, ax on Hl is given by 



where the notation =l means equality upon restriction to TLl and (Jy =l ^Wxi The independence property is useful to 

allow parallel encoded manipulations on different DFSs. However, when only a single DFS is in use, requiring independence or 
even preservation of all DFSs has no advantages, and allowing for the most general Hamiltonian as in (^) may in fact increase 
the options available for implementation. For instance, an alternative choice for encoded z and x observables is 

L _ 2 

L ^ 12' (8) 
and again cr^ =l i[(T^, crf']/2. 

In principle, based on standard universality results [|l]l, it is possible to generate any encoded unitary transformation by 
appropriately alternating evolutions under two Hamiltonians with the correct symmetry. For instance, this is certainly true if 
one can turn on/off a pair of Hamiltonians in A'^ such as, say, al and the exchange interaction E12 — (c^ • cr^ + l)/2 - for 
i[Ei2, (jI]/2 gives an encoded ay and then —i[Ei2, c^]/2 gives a^ as in (Q). 

Once encoded single-qubit manipulations are available, then universal encoded computation over DFS qubits requires the 
additional ability of implementing a non-trivial encoded gate between two logical qubits. For instance, a controlled-rotation 

gate could be constructed from a logical phase coupling of the form az^ az^' , which is supported by the natural couplings of 
many NMR and NMR-like Hamiltonians. We focus here on the first step of this program i.e., to obtain reliable single-qubit 
DF manipulations compatible with the constraints that QIP implementations unavoidably face in terms of both the form and the 
tunability of the available control Hamiltonians. 



III. CONTROLLING ENCODED QUANTUM INFORMATION WITH REDUCED ERROR RATE 

A. Two-spin NMR QIP as a case study 

The total system Hamiltonian we consider, Hs, is the sum of a time-independent internal Hamiltonian, Hint, and a time- 
dependent external Hamiltonian, Hext- The internal Hamiltonian, composed of spin-field and spin-spin interactions, is [ pT| ] 

H^nt = T^(yx<j\ + V2^l + J ■ a'^ 12) , (9) 

where i^i, ^2, and J are the chemical shifts and the coupling constant, respectively. The external Hamiltonian, describing the 
interaction between the spins and an applied RF field has the form [ pTp9[ ] 

H,xt = J2 e-'("^-^*+'^)"J/2(a;a^/2)e'('"^-^*+'^)"J/2 , (10) 

fc=l,2 

the transmitter's angular frequency lurf, the initial phase cj), and the power lu being tunable over an appropriate parameter range. 



4 



The implementation of an arbitrary unit^y gate is accomplished by modulating Hint via an external control sequence. While 
sequences can be optimized numerically p^], average Hamiltonian theory (AHT) [ ^ ] provides a systematic method for de- 
scribing any unitary propagator U{T) resulting from the evolution under the time-varying Hamiltonian Hs — Hint + H^xt in 
terms of an effective Hamiltonian H applied over the same time interval: 



T _ 

-iHT 



U(T) = Tcxp (~ i J dTHs{T)j = e 



where T is, as usual, the Dyson time-ordering symbol. In particular, AHT underlies the design of coherent refocusing and 
decoupling methods, which are able to effectively turn on/off selected contributions to the average propagator over some time 
interval. These methods have been recently revisited within the QIP context in |^ 23, n|]. We recall that the basic idea is to 



subject the system to a cyclic train of pulses V — {Pj}jLi, H^L^Pj — 1 which, in the simplest setting, are assumed to be 
infinitely short and equally spaced by > 0. The net controlled evolution over the period T = MAt can then be expressed as 

_ M 

fc=0 

where the "toggling-frame" Hamiltonians H^ are determined as Hk ~ U^HintUk, in terms of the composite pulses Uk — 



n*^^jPj , fc = 1, . . . , M, C/q = 1 In the limit of sufficiently rapid control, H simply approaches 

H=j7Y.Hk = jjY.UlH.-tU,. 



M ^ M 

k=0 k=0 



By appropriately designing the pulse sequence V, undesired contributions to H can be effectively turned off. For instance, a 
train Vi of equally spaced, simultaneous tt^ pulses on both spins (tt^tt^ pulses) averages out any phase evolution due to the cr^ 
terms in (||). Similarly, the cr^cr^ + fyCT^ coupling can also be averaged to zero by a pulse sequence V2 consisting of repeated, 
equally spaced Tr^Tr^ pulses. 



B. Universal gates via encoded dynamical control 

While AHT represents a powerful tool for designing logic gates over physical, un-encoded degrees of freedom, a direct 
application on DF encoded qubits does not automatically result in DF manipulations. Even though ensuring that H has the 
general form (||) (for instance, H £ A'^) leaves the system in a DF state, there is no guarantee that the control path has remained 
within the DFS at all intermediate times - possibly re-introducing exposure to noise. 

We begin by noting that Hint can be rewritten as 

Hint = ^{^^^^J^ + ^^W. + {^l - ^2)0': + Jd^) , (11) 

in terms, for instance, of the encoded observables - which makes it explicit that Hint G -A!^. Since both and o\a1 are 
constant on the code subspace TLl, they can be ignored and Hint further simplifies to 

=i 7r(Ai/fTf + Ja^) , Av = ~ ■ (12) 

Thus, the natural evolution implements a non-trivial logical operation within TLl. The challenge is to extract the required 
controlled operations by remaining, ideally, always within the DFS. 

The situation is simpler in the limit where, as above, control pulses are treated as instantaneous. Because Hint G one 
can ensure that each toggling-frame Hamiltonian Hk also remains in A'^ by choosing pulses such that either [f//c, Jz\ = or 
{C/fc, Jz} ~ 0. The latter condition is satisfied, for instance, by the above-mentioned pulse sequence V2, which thus implements 
a net encoded identity in this idealized scenario. 

Of course, the duration of real-life pulses is necessarily finite, and one needs to pay additional care to what happens during 
the pulse length [|T]]. In principle, DF logical operations can still be effected if sufficient control over the parameters Ai^, J in 
( p^ ) is available. The general idea is to concatenate AHT with the underlying DF encoding i.e., to implement refocusing directly 
with encoded rotations. Let us look at our DFS qubit (a more expanded account will be provided elsewhere; see also [ |3^ ] for 
related work), and imagine that encoded tt^ pulses are available as tt^ y ~ exp{~i7Ta^ y/2). Then a sequence of equally spaced 
encoded tt^ pulses (in this case a Carr-Purcell sequence [pO|]) can be used to refocus the encoded phase evolution and only leave 
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the encoded coupling active in (p^. This can be thought of as a logical or encoded "spin echo" [|3]]. A similar procedure 
holds for extracting the encoded cr^ Hamiltonian if encoded tt^ pulses are employed instead. Thus, the same schemes that are 
effective at turning on/off unwanted terms in the physical qubit evolution are effective at turning on/off unwanted terms in the 
encoded qubit evolution, provided ordinary control pulses are replaced with encoded ones. More generally, a group-theoretical 
framework extending the un-encoded approach of |p2p to encoded dynamical decoupling can be constructed. For our system, 
this implies that the ability to apply a single Hamiltonian with the correct symmetry (e.g., <t^) suffices, in principle, for gaining 
universal control. 

Unfortunately, such control is not directly available in practice, as the evolutions induced by the external RF Hamiltonian 
( p^ ) do not resemble, in general, evolutions under logical Hamiltonians. The approach we take results from the following 
compromise: we mimic the implementation of a fully encoded refocusing scheme by using available pulses whose propagator 
(not Hamiltonian) equals the required encoded rotation; we then compensate for the residual exposure to noise by control design. 
If pulse durations are optimized, then the system will reside in the DFS for a dominant portion of the computational time. In 
addition, pulse design can add robustness against noise |p4|], reducing its impact while the system resides outside the protected 
space. While in the limit of weak noise with arbitrarily long correlation times these techniques provide robustness, for realistic 
noise models actual improvements will depend heavily on the noise parameters. An explicit implementation will be reported. 

As already noted, these ideas open the way for manipulating more than a single encoded qubit. If, for instance, two DFS 
qubits are supported by the zero-quantum subspaces of, say, two proton and two carbon spins, the overall internal Hamiltonian 
will be expressible, to high accuracy and for a wide class of spin-spin coupling distributions, in terms of both single-qubit 
encoded observables ax'''^, Oz^'^ and the two-qubit encoded interaction a^'^a^^ . Thus, the ability of separately controlling 
each encoded qubit via encoded refocusing, combined with the presence of the logical phase coupling, implies the potential of 
effecting universal quantum logic with reduced error rate. 



IV. EXPERIMENTAL OUTLINE 

Liquid state NMR QIP techniques have been extensively discussed in the literature and only the salient points are recalled 
here. Because the system exist in highly mixed, separable states, NMR QIP relies on "pseudo-pure" (p.p.) states whose traceless, 
or deviation, component is proportional to that of the corresponding pure state. The identity component of the density matrix is 
unobservable and is treated as a constant under the assumption of unital dynamics {i.e., dynamics that preserves the completely 
mixed state). In this case, the evolution of a p.p. state is equivalent to the corresponding pure-state evolution. Initialization 
of the two-spin system into an intended p.p. state was accomplished using gradient-pulse techniques as described in [21 35 1 



Throughout the experimental implementation, all deviation components were explicitly verified by state tomography [36|. A 
fixed amount of identity component that optimizes the fidelity between the experimentally determined and a desired reference 
p.p. state, 1 00) (00 1 , was added to each reconstructed deviation density matrix. For each experiment, we prepared one of the p.p. 
input states = IV'miO) (i/'inOl, with y^in) {ipm \ providing a complete set of one-bit density matrices so as to allow quantum 
process tomography reconstruction [ pTHlq ]. 

Our physical system is an ensemble of Dibromothiophene molecules (Fig. 1) in a solution of CDCI3. Measured values for the 
relevant parameters are listed in the caption of Fig. 1. The experimental procedure begins with the data qubit 1 containing the 
state \ipin) to be protected, \tpin) = cq\0) + ci |1), and the ancilla qubit 2 initialized to |0). Encoding of the initial input state to 
the code space TYl is accomplished by the unitary transformation 

Uenc{co\0) + Cl|l))l|0)2 = Co|Oi) + Ci|U) , 

where Uenc is a controlled a.^ rotation on bit 2 conditioned on bit 1 having the state |0). Next, an encoded operation is performed 
on the system in the presence of noise. The information is retrieved by applying a decoding transformation Udec = C^lnc' 
producing a general output state of the form 

— Utarget |^»)l(V'm|C/tUe*® 10)2(01, (13) 

for a target single-qubit unitary transformation U target on the data spin. U target — 1 corresponds to storage of the quantum data 
under either engineered collective dephasing or natural noise, while Utarget is a non-trivial desired rotation for demonstrating 
universal quantum logic. All experiments were carried out on a 400 MHz Bruker AVANCE spectrometer 



A. Unitary and non-unitary control 



The un-encoded gate operations involved in the encoding and decoding networks were mapped into ideal pulse sequences 
usin g st andard methods [|T]|. Pulses were then implemented by modulating Hint with external RF fields as mentioned in Sect. 
Ilia pi]. 
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Non-unitary evolution, either for p.p. state preparation or for emulating collective decoherence, were implemented using 
pulsed magnetic field gradients. Magnetic field gradients take advantage of the spatial extent of the sample to induce an inco- 
herent evolution. Applying a gradient WzB = dBz/dz along the axis of the static field causes a linear variation of the Larmor 
precession frequency given by the spatially dependent Hamiltonian 

Hgrad = JzJ;,VzB/2 , 

7 being the gyro-magnetic ratio of the given nuclear species. This causes each quantum coherence pke (k 7^ £) to be multiplied 
by a spatially dependent phase factor, cxp{—ijzmke'^BzS/2), where niki is the coherence order (defined earlier) and S the 
duration of the gradient pulse. In other words, each part of the sample experiences a different coherent phase error. Tracing over 
the spatial degrees of freedom, as is done in an ensemble measurement, causes this incoherent evolution to become irreversible 
when considering the spin degrees of freedom alone. While the effects of this evolution could be immediately reversed, ran- 
dom molecular diffusion causes an irreversible spatial displacement that increases with both time and the molecular diffusion 
coefficient. Applying an inverse gradient after a time delay A (diffusion time) thus results in an exponential decay of non-zero 
coherences, exp(— A/r), with an effective noise strength given by 

-^D{j\/,BmkeSf , 

T 

D being the diffusion coefficient of the sample. Note the scaling of this decoherence rate with the square of the coherence order, 
as anticipated in the derivation of Eqs. (^. 

Using these gradient-diffusion techniques, variable strength noise can be obtained by either changing the gradient strength 
or the diffusion time. It should be noted that, in both the incoherent and decoherent case, the induced phase error is collec- 
tive to an extremely good extent, deviations from a collective action being determined by the product of the gradient strength 
(approximately 60 Gauss/cm) and the spatial displacement between the two hydrogen spins (on the order of angstroms). 

B. Reliability measures for control 

As a reliability measure quantifying the accuracy of implementing a target unitary transformation (7 on a system S we invoke 
a variant of the entanglement fidelity as introduced by Schumacher [|8||. In Schumacher's notation, let R be an auxiliary 
"reference" system, and let the initial entangled state [vj/^-^^ of the pair RS be subjected to the overall evolution 1^ (g) £^ . 
Starting from p^-^ = I'J'-'^'^) (^P^-^l, this produces a final state p^-^' = (1^ (g) (^'^'^1). Then the entanglement fidelity 

of the process £^ relative to the initial state of 5* alone, — TiR{\'i'^^){'i>^^\}, is defined as 

i^e(p^,f^) = Tr||«'«^)(*^'5|p^^'} =Tr{|*^'5)(^'«^| (l^ ® , (14) 

i.e., Ff, measures the fidelity between the input and output states of the joint system: F^, = Tx{p^^ p^^ }. can be expressed 
in terms of quantities intrinsic to the system alone once an operator-sum representation for £^ is available. If £^{p^) = 
^'Ip^^^ ^ Schumacher showed that [|3|] 

F,(p^f^) = 5:|Tr{p«A^}|\ (15) 

Because F(^{p^ ,£^) = 1 if and only if p^^ = [|3^], F^ naturally quantifies the preservation of quantum 

information - perfect preservation corresponding to implementing £^ — \^ . In the presence of the target transformation 
U = ^ 1^ , the appropriate measure should equal 1 if and only if p^^' = \J^\^^^)(^^^\U^'^ . Thus, (|l^ is generaHzed 
to a gate entanglement fidelity as follows: 

FeiU^p'^U^lS^) = Tr|c/'^|*^^)(^'^^|C/^^ (1^ (g) . 

By using the above operator-sum representation for £^ {p^), one can derive the equivalent expression 

Fe (f/ V f/'' ^ f ) = Tr { I vl/^^^) (*«^^ I ( 1 « ® £^ ) ( I vI/«^^) (*^-^ I) } = Fe (p^ , , 

where the modified dynamical map £^ — U^^S^U^ is defined by the set of transformed Kraus operators {U^^A^}. Thus, 
a perfect implementation of the desired gate corresponds to perfect preservation of quantum information under the 
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meaning of this is simply that, in the ideal case, the intended effect would be 8^{p^) = p^U^\ which is equivalent to 
ensurmg 8^{p^) = p^. Similai-to (|l5|), we then have 

F^{p^E') = Y,\lr{p'u'^Al}f , (16) 

where as above p^ = Tr/fj^'^"^) is the initial density matrix of the system alone. Taking as the standard reference state a 

maximally entangled purification |\I>^'^) for which p^ is the fully mixed state i.e., p^ — /N for a A^-dimensional state space, 
( |r^ ) finally becomes 

This form makes it explicit that the gate fidelity defined in [ ^ ] is identical with the gate entanglement fidelity formally introduced 
here. We shall still refer to the quantity in ( p^ simply as entanglement fidelity in the following. 

The above reliability measure can be related to experimentally available data, and the results take particularly simple expres- 
sions in the case of single-qubit transformations we are concerned with. Starting from the standard maximally entangled Bell 
state for the joint RS system, where S and R are now two qubits, and assuming that the process £^ = £ actually implementing 
U is unital and trace-preserving, one finds 

Fe^^ {Fu\o) + Fu\+) + Fu\+^) - l) , (18) 

where -Fc/i^^^) = Tr{J7|'!/'in)('0ir! |C/^'^(IV'Mi)('0in|)} for a generic one-bit pure input |-0i„), and |0), |+), | + i) are eigenstates 
with positive eigenvalue of az,ax, Cy, respectively. This expression was used in [ p^ for the special case U = 1. As a further 
remark, it is worth noting that as given in (|l8|) is related to the so-called average gate fidelity F proposed in via 
F = 2/3Fe + 1/3. Eq. ( p^ is directly applicable to quantifying the accuracy of DF unitary manipulations as given, upon 
decoding, by (p^). 

V. DEMONSTRATION OF A DECOHERENCE-FREE QUBIT 

The utility of a DFS to preserve quantum information (i.e., to implement the identity operation) is demonstrated under the 
action of different classes of both engineered and natural noise: a variable-strength engineered decoherent noise, a full-strength 
(crusher) engineered incoherent noise, and the natural ambient noise due to relaxation. A significant improvement in the en- 
tanglement fidelity for each class of noise is seen. In addition, because the measured entanglement fidelities remain above the 
threshold value 0.50 [p^, all implementations guarantee, in principle, the ability to preserve entanglement over a wide range of 
noise strengths. Unlike the case of an NS, the entire state of the system (data plus ancilla) remains unchanged under the action of 
the noise. While this was experimentally confirmed to a good accuracy [|4l|], we report the preservation of quantum information 
between the desired input state and the measured output state of the data qubit alone. 

A. Engineered noise 

Gradient-diffusion techniques were used to implement variable-strength noise. In order to isolate the effects of the applied 
noise, the time delay between encoding and decoding was kept fixed. In the first half of this time delay collective phase noise 
was applied to the system. Unwanted evolution due to the internal Hamiltonian was refoucsed during the second half of the delay 
by a pair of tt pulses V2 given in Sect. Ilia. The gradient strength was varied over the full dynamic range of the spectrometer 
(0 to 60 Gauss/cm) and the diffusion time A was set in such a way that a significant amount of information was lost when the 
un-encoded data spin was directly exposed to noise. This was obtained by running separate experiments with encoding/decoding 
sequences turned off. 

The experimental data are collected in Fig. 2. For the encoded data, assuming no additional loss of information with increasing 
noise strengths is consistent with the experimental results. Fitting the data with a constant value yields a value of 0.97 ± 0.01. 
The deviation from unity is consistent with the observed F^. value for the reference situation of no applied noise (i.e., the one 
implementing a net identity evolution between encoding and decoding). Therefore, these losses are caused by imperfections in 
the applied pulses as well as by natural noise processes whose action is not in the correctable error algebra Az (see below). For 
all but the smallest noise values, a substantially increased amount of quantum information is retained using the DFS memory 
than leaving the system un-protected. 
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To further confirm the robustness of the DFS memory against arbitrary noise strengths, an incoherent implementation of all 
possible collective phase errors was also realized to emulate noise in the strong dephasing limit. A single magnetic field gradient 
pulse with maximum strength was applied, causing spins on the fringe of the sample to evolve through more than 750 cycles and 
therefore acquire large phase errors. Again, the loss of information in the presence of this crusher noise is compatible with the 
measured loss due to just encoding and decoding (see Table 1). 



B. Natural noise 



The behavior of both the DFS-encoded and the un-encoded data under ambient noise was also probed in a separate series of 
experiments, with the goal of gaining qualitative insight on the relevance of fully correlated dephasing in the naturally occurring 
phase relaxation processes. In this case, the holding time between encoding and decoding was varied to allow for a variable 
exposure to noise. Because natural relaxation takes significant contributions from Ti processes (that are both amplitude and 
phase damping) in addition to transverse T2 relaxation, the unitality assumption invoked in deriving the expression (|l8|) for the 
entanglement fidelity is no longer accurate. While Fe could be still evaluated directly from ( [l7| ) upon experimentally extracting 
a set of Kraus operators, a simpler coherence metric is appropriate if the actual amplitude decay is of no interest. Similar to JT^], 
the amount of quantum coherence C (phase information) that is retained in the course of the noisy evolution can be quantified 
by experimentally determining the average off-diagonal component present in the output density matrix. Thus, corresponding to 
the two experimentally prepared transverse p.p. states |+) , |i) defined above, we calculate 

C = \{Tr{a.E{\+){+\)} + T,{ayE{\i)m) , 

where the map £ now corresponds to the natural noisy dynamics. 

The experimental data are presented in Fig. 3. Holding times ranging from a fraction of a second up to a time scale comparable 
to T2 were explored. An appreciable decay of the DFS qubit is seen in this case, witnessing the presence of non-collective phase- 
damping processes in the ambient noise. In spite of this non-robust behavior, the DFS is still able to retain quantum coherence 
much longer than the un-encoded state. This implies that a significant contribution to the overall phase relaxation is actually 
caused by fully correlated dephasing, consistent with the physical intuition based on the geometrical and chemical structure of 
the molecule. 



VI. DEMONSTRATION OF UNIVERSAL CONTROL OVER A DECOHERENCE-FREE QUBIT 



As stated in section lie, evolution according to two non-commuting encoded Hamiltonians is required for universal control. 
In the implementation, we found it convenient to adopt a choice of encoded observables intermediate between (Q), i.e. 

= —u1, cr^ — {o'x'^x + o-,V,^)/2 henceforth. In terms of this choice, and using the actual implementation parameters (see 
Fig. 1), the internal Hamiltonian (lO) is given by 

H^nt =L -137.57raf^ + S.Ttto-^ . (19) 



A. Encoded z and x rotations 

According to the above expression, the strong Hamiltonian dominates the contribution. Therefore, an encoded z 
evolution can be implemented, to high accuracy, by simply waiting the appropriate amount of time. Because of the negative 
sign in (|l9|), a positive rotation of 9 can be implemented by rotating hy 2tt — 6 about the negative logical z axis. A 
representative 7r/2 encoded rotation, cxp{—i7Ta^ /A), was experimentally implemented with a fidelity of entanglement of 0.94 ± 
0.03. Comparing this result with the accuracy of the identity operation (see Table 1), we see no significant loss of information 
due to the 2; gate. 

Implementing a x rotation is less straight forward. While a logical x Hamiltonian is present in (|l9|), it is quickly averaged 
out by the stronger logical z term. With the encoded system expressed in this form, it is clear that we must average out the 
term at a rate faster than its strength [^4|], using control operations that commute with a^. The external RF Hamiltonian 
provides this extra control parameter, at the expenses of forcing the system to leave the DFS - if only for short periods of time. 
While continuous irradiation of the spins would suffice to preserve the Hamiltonian by matching the so-called Hartman-Hahn 
condition [ ^^p3| ], as discussed in Sect. Illb a train of vr^ pulses achieves the same goal 33,4^, with the advantage that the 



system may be left in the DFS for significant portions of the control time. The trick is to note that, although the evolution 
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induced by the applied RF field does not resemble in general a Hamiltonian, for the special case of a hard tt pulse the 
resulting propagator is 

Uh = exp [~ i^icr^ + CT^) j = exp [- i-cr^ j exp [- i-a^ j = -ct^ct^ , 

T 12 r 

which mimics a net cr^ operation: the action of ?7^ ' on the code subspace is identical to the action of cr^ as a unitary operator 
(not as a Hamiltonian - note that + does not clearly respect the form (^)). 

Composite pulses |^], which provide an excellent balance between speed and robustness, were used to implement each hard 
TT pulse. Six-period pulses optimized to be robust against variations in both chemical shift (phase errors) and RF strength (control 
errors) [ p4| ] were used to emulate the required sequence of encoded tt^. Each of these hard pulses is 62.4/is in duration and is 
followed by a delay of 630/is. Therefore, the system resides in the protected space for over 90% of the computational time. The 
phases of the tt^ pulses were alternated systematically as described by a WALTZ sequence [p^], so as to minimize the impact 
of experimental errors. In particular, a 64-cycle sequence was used to achieve a 7r/2 encoded rotation, exp(— z7ro'^/4), with a 
fidelity of entanglement of 0.94 ± 0.03. Again, we see no significant loss of information due to the x operation. 



B. Composite encoded y rotation under collective phase noise 

To explicitly test the robustness of the available logical x and z manipulations, a composite encoded rotation by 7r/2 about y 
was implemented in the presence of variable-strength collective phase noise i.e., the sequence of encoded rotations 

exp (- ij'^z) exp (- ijf^x) exp {^j'^i) = exp (- ijf^y) 

was performed. As in the DFS memory experiments, gradients were used to induce a spatially incoherent error over different 
parts of the sample. However, the noise effects associated with a time-independent gradient Hamiltonian (i.e., an infinitely long 
correlation time) would tend to be effectively averaged out over time by the applied control sequences as described by coherent 
averaging [|o|| and dynamical decoupling [|2|] . In order to make sure that the net action of the applied gradients is maintained in 
the presence of the external control, a procedure similar to the fast-switching control schemes discussed in [ ^3| ] was followed, by 
rapidly modulating the strength of the applied gradient Hamiltonians over the course of the control sequence. Thus, a temporal 
incoherence was also superimposed at each spatial location in the sample, enforcing a finite correlation time Tc (hence a non-zero 
cut-off frequency) in the spectral density describing the noise. In practice, the gradient waveform was determined via a random 
walk process, whose shape is depicted in Fig. 4. The gradient strength was changed every 50.6 /is therefore Tc ^ 50.6/is. By 
making sure that Tc is short compared to the control cycle time of the sequences used to implement the composite rotations 
(~ 700 /is in our case), active control is made ineffective at averaging out the high frequency effects of the noise during the 
computation. 

A broad range of values for the maximum applied gradient strength were explored to test the robustness of the computation 
to collective phase errors. The experimentally determined gate entanglement fidelities are shown in Fig. 5. As expected, the 
computation is protected up to a particularly noise intensity and then falls off with increasing noise strength. As in the memory 
case, it is worth stressing that values well exceeding the value 0.50 have been achieved over the entire range of applied noise 
strengths. It should be noted that because the active sample is order 1 cm, most of the sample is experiencing noise strengths 
significantly stronger than natural fluctuations - which are approximately 1 Hz in strength. 



VII. CONCLUSIONS 



We have provided the first demonstration of universal control over a DFS-encoded qubit. The implementation relied on 
combining the benefits of passive noise protection via DFS coding with the ability of relaxing the constraints of fully DF manip- 
ulations via appropriate control design - thereby also validating the underlying principles of encoded decoupling. We believe 
that our techniques are applicable to a wide class of quantum information devices, where collective dephasing mechanisms play 
a dominant role and where the structure of the system's internal Hamiltonian can be mapped onto a NMR-type Hamiltonian. 
These may include various solid-state proposals as discussed in Thus, our results improve the prospects that DFS/NS 

coding, combined with encoded dynamical decoupling and robust control design, will play a practical role for both protected 
storage and manipulation of quantum information in QIR 
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TABLE 1. Experimental data for the implementation of full-strength collective dephasing. Input-output fidelities and entanglement fidelities 
corresponding to the application of the intended error model to both the DFS encoded (Qz,df) and the un-encoded data spin (Qz,un) are listed, 
along with the values relative to the reference situation of zero applied noise between DFS encoding and decoding (Qo,df)- Crusher gradient 
fields with full strength ^ 60 Gauss/cm were applied for a period S — 745/is. The measured values for the un-encoded test data confirm the 
expectation that the applied noise process induces full phase damping on the data spin, with predicted Fe — 0.50. Systematic uncertainties are 
~ 0.02 while statistical uncertainties are ~ 2%, both due to errors in the tomographic density matrix reconstruction. 
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FIG. 1. Molecular structure of Dibromothiophene. The two proton qubits are indicated. As spectroscopically the two protons are effectively 
undistinguishable, the qubit labels have a purely formal meaning. All experiments were carried out in a magnetic field of ~ 9.7 T with one 
proton on resonance. The frequency shifts of the second proton is V2 — 137.5 Hz, while the J-coupling constant is J = 5.7 Hz. The 
longitudinal and transverse relaxation times are Ti ~ 7 s and T2 ~ 3.5 s, respectively. 




Noise strength 1/t (s-^) 



FIG. 2. Experimentally determined entanglement fidelity for the implementation of variable-strength collective dephasing. Both the 
behavior of the DFS-encoded (squares) and the un-encoded (circles) data is shown. The independent axis (noise strength) was determined by 
fitting the un-encoded data to an exponential decay of the form Fe = A exp(— teu/r) + 0.5, with tev = A + 25 — 37.765 ms. The un-encoded 
data is only displayed for reference. The encoded data is fit to a constant value Fe — C, yielding the best estimate C — 0.97±0.01. Systematic 
uncertainties (not included in the figure) are ~ 0.02. 
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FIG. 3. Experimental data for the phase information retained after exposure to the natural system noise. The average preservation of and 
ay was measured as a function of holding times from to ~ 3 s. Improvement over the un-encoded case is seen, confirming that collective 
phase errors are one of the dominating modes of natural noise for this system. 
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FIG. 4. The temporal variation of the gradient waveform used to test control of a DFS qubit in the presences of noise. The shape was 
determined by a random walk algorithm. The stepping time was 50.6/^s, which is faster than any of the control time scales relevant in 
implementing the encoded transformations. 
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FIG. 5. Experimentally determined entanglement fidelity for the implementation of a composite encoded y rotation of 7r/2 in the presence 
of noise. Magnetic field gradients implement a spatially incoherent collective phase error as a function of the molecular position. Gradient 
strengths from to ~ 100 kHz/cm were applied over a 1 cm sample. The behavior of Fe is flat over a broad range of noise strengths and 
remains significantly above the 0.5 threshold for all noise values considered. This convincingly demonstrates the ability to control a DFS qubit 
in the presence of noise significantly stronger than the natural noise of the system. 
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